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Using the functional integral method, we construct a theory of heterotic SIS Josephson junctions 
between single- and two-gap superconductors. The theory predicts the presence of in-phase and 
out-of-phase collective oscillation modes of superconducting phases. The former corresponds to the 
Josephson plasma mode whose frequency is drastically reduced for ± s-wave symmetry, and the latter 
is a counterpart of Leggett's mode in Josephson junctions. We also reveal that the critical current 
and the Fraunhofer pattern strongly depend on the symmetry type of the two-gap superconductor. 
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The Josephson effect is one of the most drastic phe- 
nomena in superconductivity [1]. Cooper pairs can tun- 
nel through an insulating barrier in a non-dissipative 
manner. This particular feature has attracted tremen- 
dous attention of not only physicists but also device en- 
gineers. 

Very recently, multi-gap superconductors have been re- 
visited since the discovery of an iron-based high-T c su- 
perconductor @, 0, H|. In contrast to cuprate high-T c 
superconductors, 3-d electrons on the iron atom form 
multi-bands whose Cooper pairs condense into a multi- 
gap superconducting state. The angle resolved photoe- 
mission spectroscopy has reported that each of multiple 
disconnected Fermi surfaces is fully gapped [5] and other 
experiments have also supported the gapful features [6[. 
On the contrary, the nuclear magnetic resonance have 
shown typical gapless features [7J. 

In order to compromise the controversy, the presence 
of ± s-wave gaps on the disconnected Fermi surfaces has 
been proposed @, 0, [To[. The essence of the ± s-wave 
symmetry is a sign change between different s-wave order 
parameters. This is expected to bring about novel behav- 
iors in phase interference effects. In particular, Joseph- 
son effects in SIS junction between the single- and the 
± s-wave multi-gap superconductors as schematically il- 
lustrated in Fig. [U drastically reflect the sign change. In 
this Letter, we focus on such a heterotic SIS junction 
and clarify peculiar Josephson effects. We have three 
main results, i.e., the drastic reduction of (i) the Joseph- 
son plasma frequency, (ii) the critical current, and (iii) 
the Fraunhofer pattern visibility. The ± s-wave symme- 
try leads to a cancellation between the two Josephson 
currents which arise from the two tunneling channels in 
this system. 

In the proposed junction as shown in Fig. [TJ the left 
(right) electrode is a single- (two-) gap superconductor 
with the width sl (sr). The insulator width and the 



dielectric constant are d and e, respectively. The cur- 
rent and the magnetic field are applied along z- and y- 
direction, respectively. Similar situations were also ex- 
amined from other viewpoints [111 ]. 



The system's Hamiltonian H — f R d 3 r(H^ + H 



^pair^ _j_ jj^^ w here TC R describes 

the kinetic energy of the i-th band electrons in the 
right electrode. The pairing term in the right electrode 
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^i2(^ 1)+ ^ 1)t ^f +h.c.) [12], where g t > 0, and the 
inter-band interaction can be either attractive (i.e., gi2 > 
0) or repulsive (i.e., gi2 < 0). The left electrode is de- 
scribed by Ti^-^H^ 11 , where Hf^ is the kinetic energy and 
^pair = — g s ip^ ip s ^ Tp*^, in which g s > 0. The tunneling 

Hamiltonian Ht = + where means the 

tunneling between the electrons in the left side and the 
i-th band electrons in the right side. Using the imaginary 
time functional integral method [13, LL4[ , the effective ac- 
tion with respect to the order parameters AW and A s is 
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FIG. 1: A schematic figure of the present heterotic junction 
system. The left electrode is a single-gap superconductor, and 
the right electrode a two-gap superconductor. 
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given by S eff = J^dr[f R d 3 r(^| A« | 2 /# + 9i I A® | 2 /# + 
^n)+/ L <i 3 r|A s | 2 /^ s ] -Tr In G -Tr In G -1 , where /? is the 
inverse temerature and T^ n = — #12 (A^ 1 )* A^ 2 ) + c.c.)/g. 
We assume here that g = gig 2 — g\ 2 > 03 • The Green 
functions for the non-interacting system Go and the to- 
tal system G are 6x6 matrices. We do not write their 
explicite expressions here, but those consists of 4 x 4 for 
the two-gap and 2 x 2 for the single-gap superconduc- 
tors [13, ll4|. The inter-band Josephson coupling term 
V [n is rewritten as V [n = -2(g 12 /g)\A^\ |A^| cos^ 1 ) - 
(pW), in which AW = |A«|e^ (i) . 

Based on the standard procedure the effective 

Lagrangian density of the superconducting phases on zx- 
plane in the real time formalisum is given by 



with respect to A? and Af as follows, 
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where the dimensionless parameters in each electrode are 
defined as a = e/i 2 / s^d, cti = e/i 2 /sR<i, a~ x = a^ 1 + 
a 2 x , r] = A 2 / ' s\,d, rji = A^) 2 /sR<i, and f/ -1 = 77^ 1 
The magnitude of the electric (magnetic) field coupling is 
characterized by a and ai (77 and rji) [l5[. The constants 
C and L are defined as C = 1 + a + o; and L = 1 + 77 + 77, 
respectively. Equations (|6]) and Q correspond to the 
generalized Josephson relations [16| . The Euler-Lagrange 
equation with respect to A RL gives the Maxwell equation, 
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and note that ao = (h/e*)dtip s + A^, a x = (hc/e*)d x (p s — 
At, a ( = (fi/e*)ft^ (i) +A0 , a, (i) = (fic/e*)^^)-^, 
and e* = 2e. The phase <p s is defined as A s = |A s |e^ s , 
and ji is the Josephson critical current between z-th and 
single-band Cooper pairs. The charge screening length 
and the penetration depth on the left (right) electrode are 
\± (/i^) and A (A^), respectively. The last term in the 
gauge-invariant phase difference (|4|) is the z component 
of the spatial averaged vector potential in the insulator, 



defined as Al 



dr 



*rl — u 1 f-d/2 A z (z) dz. The electric and 
the magnetic fields in the insulator are defined as E RL = 
-c~ l d t A RL - cT 1 « " O and B« RL = d -\A R -At)- 
d x A RL , respectively. Here, let us focus on the Josephson 
coupling energy ([2]). The first and the second terms are 
the ordinary Josephson coupling terms, while the third 
term corresponds to the inter-band Josephson coupling 
energy and | J[ n \ is proportional to |(<7i2/#)sr|. One finds 
that Ji n is positive (negative) if gi 2 > (#12 < 0). If 
sl and sr are much larger than d, then it allows us to 
regard | J in | > j u j 2 . 

From Eq. (pp), we have the Euler-Lagrange equations 



where A 



Ji 



Air e*dji/ he 2 . The first term on the right 
hand side of Eq. (|HJ) is the summation of the Josephson 
current terms [Fig. [T]. Using Eqs. (|6|)-(|8|), we obtain 
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Next, from the Euler-Lagrange equations about cp s and 
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where the parameter £ (£) means the difference of the 
magnitude of the electric (magnetic) field coupling be- 
tween the different superconducting bands as £ = [a\ — 
0L2)/ (ai + a 2 ) and (" = (771 - 772)/ (771 + 772). 

The present description is valid from very thin elec- 
trode junctions (sl ~ M and sr ~ //W) to conventional 
thick ones. In the latter case (sl ^> M and 5r ^> /i^), we 
can take an approximate treatment, a — > and — > 0. 
Remark that the present paper highlight, i.e., particular 
features due to ± s-wave is unchanged in this limit. 

Now, let us examine the collective modes involved in 
Eqs.© and (|TQ|). For this purpose, we linearize them 
around a stable point of V. First, we focus on J- in > 
0. Then, a stable point for (6^,6^) is (0,0) and the 
dispersion relations uS±(kx) is given by 
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where 



<4=C(u>l 1 +u>l 2 )[l + L-\k x /Kn 
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The Josephson plasma frequency associated with the 



Josephson current for 



c/y^Aji, while the 



pseudo Josephson-plasma frequency associated with the 
inter-band Josephson current, u- m = ^/ Aire* d\J- in \ / eh. 
Note that K 2 = c _2 e(cjp 1 + cj 2 2 ), the dimensionless pa- 
rameters X and Y are defined as X = 1 — £ 2 C _1 (1 + 
a) and Y" = 1 — C 2 L _1 (1 + 77), respectively, and the 
quantities D and R are, respectively, D = C^oA + 
cu^i-S+CL-Hkcc/K) 2 ] and i? = X{(l-£ 2 )( w 2 - w £) 2 + 
AaC-^D-^+uji)^} 2 }, where 5 = (ji -j 2 )/(ji+j 2 ). 
When £ = £ = 5 = (i.e., the superconducting charac- 
ters are perfectly equivalent between the two bands), we 
find that 00+ = uop and uo- = uo\,. We then notice that 
no term related to the inter-band Josephson coupling is 
involved in the expression of ujp. It indicates that the 
origin of cop is irrelevant to the motion of the relative 
phase if. Then cj + corresponds to the in-phase motion 
for 6^ and 0^ . On the other hand, the origin of uj- is 
the out-of-phase motion for 9^ and 0^ 2 \ Next, we study 
another case of J- in < 0, in which (7r, 0) is a stable point 
since |Ji n | ^> ji,j2- Expanding V around (7r, 0), we have 
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(L^ikx/K) 2 ) and i?' = X{(1 
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and (b) show the typical dispersion relations for J- in > 
and Ji n < 0, respectively. For both cases, the frequency 
of the out-of-phase mode co- is found to be lower than 
the in-phase mode o;+ for an arbitrary value of k. Here, 
we take the limit k x — > in Eqs. ([TT]) and ([T2j) to explic- 
itly evaluate the gap frequency for these modes. As for 
Jin > 0, the leading order terms are, respectively, given 
as u+ ~ (cjp! + oo 2 2 ) 1/2 and cj_ ~ {a x + a2) 1/2 ^in by 
regarding a and a to be small. However, remark that 
we keep the term au^co^ + ^ 2 2 ) _1 in the above eval- 
uation, because | Ji n | > ji, 32 even though a is small. 
Similarly, when Ji n < 0, we have o;+ ~ la; 2 -, — a; 2 -! 1 / 2 



and cj_ 



(ai + a 2 ) 1/2 ^ 



The gap of o; + is char- 



acterized by a superposition (Ji n > 0) or subtraction 
(Jin < 0) between oo p i and cu P 2. Thus, we refer cj + to the 
Josephson plasma mode. We emphasize that the a sig- 
nature of ± s-wave is the reduction of plasma frequency 
[Fig. [2]. On the other hand, since the gap of ou- is char- 
acterized by v- lYl and c^, we find that it corresponds to 




FIG. 2: The dispersion relations for uj±. The solid (dotted) 
line is for cj+ (<^-). The parameters are set as follows: a — 
OLi - 1CT 3 , a 2 /ai - 1.4, rj = rji = 10 3 , 772/771 - 1.4, j 2 /ji - 
0.9 and \J- in \/ji — 40.0. The reduction of u;+ is observed for 

Jin < 0. (a) Jin > 0. (b) Jin < 0. 



the gap of the Leggett's mode, which was derived as a 
collective mode generated by the density fluctuation be- 
tween two superfluidities [17]. Thus, it should be called 
Josephson-Leggett mode. The inter-band Josephson cou- 
pling and the charge density fluctuation create the the 
mode. Conventionally, the first and the second terms in 
Eq. (pQ) are fixed to be zero (ctj — > 0) because the charge 
screening length is much smaller than the electrode size. 
Then, the effective Lagrangian density gives the stan- 
dard Josephson relation d t 0^ = (e* J/ft)i^ L , resulting in 
d t (0^ — #( 2 )) = 0. It means that the Josephson-Leggett 
mode becomes a gapless mode. In contrast, the mode cop 
can still remain massive in a{ — » 0. The retainment of 
non-zero is responsible for the finite gap frequency of 
the Josephson-Leggett mode. The bulk Leggett's mode is 
normally embedded in the quasi-particle excitation con- 
tinuum [l8j], while the Josephson-Leggett mode is more 
clearly and easily observable because the mode lies far 
beneath the gap energy. 

The remaining part of this Letter is devoted to basic 
Josephson effects. First, let us discuss the Josephson 
critical current j c . The bias current is assumed to be 
uniformly applied without the external magnetic field. 
Namely, 0W and 0^ are assumed to be unifrom along the 
x-axis in Fig. [TJ The bias current I is added to the right 
hand side of Eq. ([9|) with the elimination of the spatial 
dependent terms, and j c can be derived by estimating 
the maximum threshold of I under keeping a stationary 
solution. The condition is given by 
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Equation ([T3j) means that / coincides with the sum of two 
Josephson currents between the electrodes, while Eq. ([T4|) 
is an internal current conservation law, which gives a 
significant constraint on the critical current. When Ji n > 
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FIG. 3: The current vs. the magnetic flux. The value of Io is I(H V = 0) for <p = 0. (a) Ji n > and <p = 0. (b) Ji n < and 
(p = 7r. No current is observed when ji = j2. (c) The demonstration of the cancellation between two Josephson currents. 



0, the preferable choice of ip is 0. Equation ([H]) implies 
(ji - J2) sin^ 1 ) = 0, because = 9^ 2 \ This is always 
satisfied if ji = j'2. Then, since 0^ can vary from to 
2?r, j c = ji + j2. If ji 7^ j2, then can deviate from 
and j c < ji + j2- On the other hand, when Ji n < 0, 
should be tt. Equation (jT4j) implies (ji -f j^sin^ 1 ) = 0. 
The only possible solution is 0^ = and 0^ = — tt, 
because j\ + 32 7^ 0. Thus, we find that the value of j c 
is drastically reduced compared to the case of J\ n > 0, 
e.g., j c = for the case of perfectly identical ± s-wave 
two-gap superconductivity. 

Next, we consider the Josephson effects in the pres- 
ence of the external magnetic field H y . We focus 
on stationary solutions, i.e., we drop the temporal 
terms of 0W. According to Eq. jTJ), 6^ 1 \x) = kx + 
#0 + (v/ r l2) ( p(x) and 6> (2) (x) = kx + 6 - (77/771 
where /c = L(e*d/hc)H y and # G [0, 2tt) is an in- 
tegral constant. The observed current is then given 
by I(H y ,6 ) = J^ f x 2 /2 [j^mO^ix) + j 2 sm0^(x)} dx. 
Hereafter, we assume that cp(x) is spatially uniform. 
Taking account of |Ji n | > ji, j 2: we should have « 
— Ji n sin(^ from Eq. (flQ|) . When J in > (i.e., ip = 0), 
the magnetic field dependece of the current is given 

by I(HV,0 ) = L x sinfl [0"i + J2)(*o/7r$)sin(7T$/$ )], 
where $ = 27rftc/e* and $ = LH y dL x . Then, I(iJf) = 
max, \I(Hy,6 )\ = L x (ji+j2)|(*oM) sin(7r$/$ )|. As 
a result, we obtain the ordinary Fraunhofer diffraction 
pattern as a function of the magnetic flux <I> [Fig. [3(a)]. 
The maximun value of I(H y ) is the sum of two Josephson 
currents as 1(0) = L x (j\ +32)- We also observe that the 
net current conventionally vanishes when <I> = k&o (ft G 
N). In contrast, as for Ji n < (i.e., ± s-wave), the current 
is given by I(H y ) = L x \j x - j 2 ||($oA$) sin(7r$/$ )|. If 
ji = J2, the Fraunhofer diffraction pattern completely 
disappears. When ji 7^ j'2, the pattern is observable 
except for <I> = ft3>o, but the maximum value becomes 
unexpectedly small [Fig. [3(b)]. The situation at ji = j 2 
is schematically displayed in Fig. [3(c). The Josephson 
currents for 6^ and 0^ cancel out each other. 

Finally, let us discuss how to experimentally confirm 
the theoretical predictions. We point out that the maxi- 
mum Josephson current can be estimated from the nor- 
mal state resistance based on Ambegaokar-Baratoff rela- 
tion [IK under an assumption J- in > 0. If the measured 
j c is significantly reduced from the one estimated above, 
then Ji n < 0, i.e., ± s-wave symmetry is concluded. 



In summary, we microscopically derived an effective 
Lagrangian density of the SIS Josephson junction be- 
tween single- and two-gap superconductors and exam- 
ined the collective modes, the critical current, and the 
Fraunhofer pattern. We found that these properties are 
considerably affected by the type of the pairing symme- 
try of the two-gap superconductor. We conclude that the 
heterotic junction is useful to identify directly a symme- 
try of two-gap superconductors. 
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